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In the present work, we theoretically study the nonlinear regime of charge transport through a 
quantum dot coupled to the source and drain reservoirs. The investigation is carried out using 
a nonequilibrium Green’s functions formalism beyond the Hartree-Fock approximation. Employed 
approximations for the relevant Green’s functions allow to trace a transition from Goulomb blockade 
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response to thermal gradient applied across the system are discussed, including experimentally 
observed thermovoltage zeros. 
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I. INTRODUCTION 

Thermoelectric properties of mesoscale and nanoscale 
systems are being intensively studies in the last two 
decades. Tailored nanostructures such as quantum dots 
(QD) and/or molecules sandwiched in between conduct¬ 
ing electrodes attract significant research interest partly 
because they are expected to be useful in manufactur¬ 
ing of highly efficient energy conversion devices. The im¬ 
proved efficiency of heat-to-electricity conversion in these 
systems originates from sharp features appearing in their 
electron transmission spectra, as it was predicted in ear¬ 
lier studies of thermoelectric properties of solids [H-Q- 
Apart from possible applications, studies of thermoelec¬ 
tric properties of nanoscale systems can provide a deeper 
insight into the nature and characteristics of electron and 
thermal transport mechanisms Recently, a new re¬ 

search field of nanoscale thermoelectronics has emerged 
and thermoelectric properties of tailored nanosystems 
have been explored both experimentally and theoretically 

[JQ- 

As known, heat-to-electric-power converters operate 
due to Seebeck effect. The latter occurs when thermal 
and electric driving forces simultaneously affect electron 
transport through the considered system. A thermovolt¬ 
age Vth appears when the temperature differential Ad 
is applied across the unbiased system provided that the 
electric current is completely suppressed. Thus the emer¬ 
gence of Vth indicates the energy conversion. Another 
quantity characterizing thermoelectric transport through 
QD, molecular junctions and other nanoscale systems of 
similar kind is thermocurrent Ith- It may be defined as 
a difference between the electron tunnel current flowing 
through a biased system in the presence of a tempera¬ 
ture differential and the current flowing at Ad = 0. As 
well as the thermovoltage, Ith is controlled by simulta¬ 
neously acting electric and thermal driving forces, and 
the combined effect of these forces depends on the bias 


voltage polarity and on the type of charge carriers in¬ 
volved in the transport process. Assuming for example 
that the left electrode of the considered junction is kept 
at higher temperature than the right one, Ith takes on 
negative/positive values when charge carriers are corre¬ 
spondingly electrons/holes. 

When Ad <C d^^fl (d^^/j being the temperatures of 
the left and right electrode, respectively) the system 
operates within the linear in temperature regime, so 
Vth = —SAT. Within this regime, the thermopower S 
describes the efficiency of energy conversion along with 
the thermoelectric figure of merit ZT. Correspondingly, 
properties of both S and ZT in nanoscale systems have 
been (and still are) intensively studied |l,sl422| . As the 
temperature differential across the system increases, the 
system may switch to nonlinear regime of operation. For 
example, a thermovoltage that nonlinearly changes with 
Ad was observed in experiments on semiconductor QDs 
and single-molecule junctions Presently, signifi¬ 

cant effort is being applied to theoretically analyze ther¬ 
moelectric properties of nanoscale energy converters op¬ 
erating beyond linear regime 

It is known that Coulomb interactions between charge 
carriers may strongly affect thermoelectric properties of 
QD/molecules leading to novel and distinct phenomena. 
Manifestations of Coulomb interactions vary depending 
on relationship between the charging energy U, the cou¬ 
pling strength r characterizing dot-leads contacts and 
the characteristic temperature d = + On). When 

the dot/molecule is weakly coupled to the leads, so that 
U significantly exceeds both r and the thermal energy 
kO {k is Boltzmann’s constant) the equilibrium density 
of states of electrons on a single-level dot (DOS) dis¬ 
plays two peaks whose separation equals U, as shown 
in the Fig. 1. As the coupling of the dot to the elec¬ 
trodes strengthens, the peaks in the DOS become lower 
and broader. A furthermost growth of r results in dis¬ 
appearance of the Coulomb blockade peaks. When the 
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FIG. 1: (Color online) The equilibrium electron density of 
states on a single-level quantum dot. The curves are plot¬ 
ted assuming that Eo = —2meV\ U = GmeV\ at kOL = 
kOn = O.lGmel^ (left panel) and kOL = kOu = 0.6meV 
(right panel). 


temperature 0 takes on values below the Kondo temper¬ 
ature 9k , a single sharp maximum emerges in their stead. 
At higher temperatures, the DOS becomes smooth and 
featureless. 

The effects of Coulomb interactions on thermoelectric 
transport throu gh Q D/molecules were intensely studied 
(See e.g. Refs. 15, la, 13, 23,133). In the present 

work we contribute to these studies by analyzing the ef¬ 
fect of Coulomb interactions on the nonlinear thermo¬ 
voltage and thermocurrent. 


I. MAIN EQUATIONS 

In the following analysis we concentrate on electron 
contribution to thermoelectric transport omitting from 
consideration the phonon part. Retaining a single orbital 
on the dot/molecule, we write the relevant Hamiltonian 
as H = H]j + Hl + Hu + Ht- Here, the first term 
describes the dot. It is taken in the standard form: 

Hd = '^ E^dlda + Ud^dad^^^d-a- ( 1 ) 

(T 

Here, are creation and annihilation operators for 

the electrons on the dot with a certain spin orientation, 
E„ = Eq is the energy of a single spin-degenerated dot 
level and U is the charging energy. The terms Hp (/3 = 
L, R) are corresponding to noninteracting electrons on 
the left/right electrode: 

Hfj = 'y ) irP<7cl.^^Cri3a ( 2 ) 

r(7 

where Cr/jo- are single-electron energies on the lead /3 
and cl^^, CrjScr are creation and annihilation operators 
for electrons on the leads. The last term: 

Ht=J2 + H.C. (3) 

r(3a 


describes tunneling effects between the dot and the elec¬ 
trodes. The factors Trpa are coupling parameters char¬ 
acterizing the coupling of the electron states on the 
dot to the leads. For a symmetrically coupled system 

'TrLa — '^rRa = 'Tra- 

Now, we employ the EOM method to compute the re¬ 
tarded Green’s function for the dot. Disregarding spin- 
flip processes, we arrive at separate expressions for the 
Green’s functions corresponding to different spin chan¬ 
nels. These expressions can be presented in the form: 


G^^(E] — _ A Eq S]q 2 I^(I (^-g-)) _ 

" ^ {E-Eo- ^oa)[E -Eo-U- S-2) + C/Si, ■ 

(4) 

Here, {n-a) are one-particle occupation numbers on the 
dot: 


{n.) = ^ j dEHn{G<{E)). (5) 

Self-energy corrections included into the expression (|3]) 
have the form: 


_ \ ' |'A/3(t| — I 

^Ocr — y = ^Ocr + 

^ E - erPa + ^'n 


( 6 ) 


rP 


rP 


1 


E - Crp-cr + ir] 



where is the Fermi distribution function for the en¬ 
ergy Crpa- and chemical potential p.p and rj is a positive 
infinitesimal parameter. The expression (j4]) was first ob¬ 
tained by Meir et al [sll. Later, the same expression was 
derived and employed in several works where transport 
properties of quantum dots and molecules were theoreti¬ 
cally studied. 

The lesser Green’s function G^{E) is related to the 
retarded and advanced Green’s functions {G'^J{E) and 
G““(A), respectively) by Keldysh equation: 


G<{E) = GI''{E)E<{E)GT{.E). (10) 


Both charge and energy transfer through a QD/molecule 
weakly coupled to the electrodes may be successfully 
studied employing NEGF within the Hartree-Fock ap¬ 
proximation.This method brings simpler expressions for 
the Green’s functions than those given by Eqs. dll)-®- 
Being employed to study current-voltage characteristics 
in a single-level quantum dot, it yields typical stair-like 
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curves with two steps maintaining the correct 2 : 1 ra¬ 
tio of heights of two successive steps for a dot symmet¬ 
rically coupled to leads [s^. However, the Hartree-Fock 
approximation cannot be used to analyze charge/energy 
transport through a QD/molecule strongly coupled to the 
leads at low temperatures. To catch the Kondo peak and 
study related transport phenomena one needs to compute 
Green’s functions beyond the Hartree-Fock approxima¬ 
tion. To correctly describe transport within the Kondo 
regime, it is important to choose a suitable approxima¬ 
tion for the lesser Green’s functions. 

In further calculations, we approximate as fol¬ 

lows: 


E<iE)=^Y,f^iE)^P{E) ( 11 ) 

/3 

where r^(£l) = —2Im[E°^(£l)] and f^{E) is the Fermi 
distribution function for the left/right electrode. This 
approximation leads to the correct 2 : 1 ratio of heights 
of two steps displayed on the current-voltage curves for 
a symmetrically coupled system (F^ = F^ = F) within 
the limit of weak coupling (F <C H), as shown in an 
earlier work 0- On these grounds, we believe that 
Eqs. (Hl)- (fTT]) may be employed to study characteristics of 
thermoelectric transport through QD/molecules in both 
Goulomb blockade (weakly coupled systems) and Kondo 
regime (strongly coupled systems below Kondo temper¬ 
ature) . 

For a symmetrically coupled system including a single¬ 
energy QD/molecule as a linker, electric current is given 
by the Landauer expression: 

^ = ( 12 ) 

Here, T{E) is the electron transmission function: 

T{E)=^-TJ2[G7{E)-GT{E)]. (13) 

(T 

where the wide-band approximation is used for the cou¬ 
pling parameter F. The Green’s functions included in 
(fOl) implicitly depend on temperatures of the electrodes 
0L and Or because self-energy terms Eio- given by 
Eq. ([7]) include temperature dependent Fermi distribu¬ 
tion functions Also, in the presence of temperature 
differential AO = 9l — Or the bridge energy Eq may 
be renormalized including a correction proportional to 
AO namely: Eq = Eq + zkAO where z is a dimen¬ 
sionless factor whose magnitude is of the order of one 
or smaller [33| . Therefore, the electron transmission de¬ 
pends of the temperatures Or, Or. In the following anal¬ 
ysis we assume for certainty that the right electrode is 
kept at a fixed temperature, and it remains cooler than 
the left one, so AO > 0. Basing on Eqs. (fT])- (fT^ one 
may theoretically analyze characteristics of thermoelec¬ 
tric transport through the considered systems. 


III. THERMOVOLTAGE 


As known, the thermovoltage is determined from the 
open-circuit condition I{Vth,A0) = 0 assuming that no 
bias voltage is applied across the system. Solving this 
equation one obtains Vth{A9). In the present work Vth 
was computed at two different values of 9r. These val¬ 
ues were chosen in such a way that one of them is below 
the Kondo temperature whereas another one is above the 
latter. Obtained results are presented in the Fig. 2. At 
small values of temperature differential the thermovolt¬ 
age is a linear function of AO which can take on positive 
or negative values depending on the nature of charge- 
carriers (electrons/holes) involved in transport. As A9 
increases, Vth reveals a distinctly nonlinear dependence 
of A9. 


With increasing AO the thermovoltage magnitude in¬ 
creases, and it reaches its maximum/minimum at a cer¬ 
tain value of AO. Further increase of the temperature dif¬ 
ferential reduces Vth magnitude. For every value of the 
bridge energy Eq, there exists the corresponding value 
of AO were Vth becomes zero. Subsequent growth of 
the temperature differential leads to emergence of ther¬ 
movoltage with the reversed polarity. We remark that 
some curves plotted in the Fig. 2 qualitatively agree with 
those obtained in recent experiments on semiconducting 

QD m. 


The reversal of Vth polarity may be explained as fol¬ 
lows [i^. Assuming for certainty that chemical poten¬ 
tials of the electrodes in the unbiased system ^r = fJ,R = 
0, and the charge carriers are electrons, the difference in 
the temperatures of electrodes leads to the flow of elec¬ 
trons from the left (hot) to the right (cool) electrode. 
To suppress this thermally excited current, a negative 
thermovoltage emerges which grows in magnitude as A9 
increases. However, as Or rises, the step in the Fermi 
distribution function for the left electrode is being par¬ 
tially smoothed out. This opens the way for holes to 
flow to the right electrode. At certain value of tempera¬ 
ture differential the hole flux completely counterbalances 
the electron flux. So, at this value of temperature differ¬ 
ential the thermally excited electric current vanishes at 
Vth = 0 . 


As known, the channels for thermally excited electron 
transport are opening up when a maximum in the DOS 
is located near the chemical potential of the electrodes. 
When we consider the system within the Kondo regime, 
the DOS reveals a single maximum near E = 0 which 
is shown in the left panel of the Fig. 1. So, the elec¬ 
tron transport occurs through E = 0 transport channel. 
However, both shape and position of the Kondo peak 
depend on the position of the bridge energy level Eq. 
Also, in the presence of a temperature differential, Eq 
becomes renormalized itself. It appears that even slight 
shifts in the bridge level position may result in significant 
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FIG. 2: (Color online) Thermovoltage as a function of the 
temperature differential for a system below (left panel) and 
above (right panel) the Kondo temperature presented at sev¬ 
eral values of the bridge energy Eq. Left panel: The curves 
are plotted at 17 = QmeV, T = 3meV, ROr = Q.lQmeV (left) 
and r = O.SmeV) ROr = O.GmeV (right). 



FIG. 3: (Color online) Left panel: The effect of coupling of 
the electrodes to the bridge on the thermovoltage. The curves 
are plotted at 17 = 6meV, Eo = —ImeV, R6r = O.GmeV. 
Right panel: The thermovoltage as a function of temperature 
differential at several values of charging energy U. Curves are 
plotted for R6r = O.GmeV; Eo = — lmeI7; F = 0.3meV. 


variations of the intensity of electron transport through 
the channel associated with the Kondo peak. The more 
suitable for transport this channel becomes, the greater 
values Vth takes. Comparing the curves displayed in the 
left panel of Fig. 2, one observes that Eq = —ImeV 
provides the most advantageous conditions for the ther¬ 
mally excited transport below the Kondo temperature 
than other values of Eq used to plot these curves. 

When the temperature of the coolest region rises above 
0k, the Kondo peak in DOS disappears but two broader 
and more robust peaks corresponding to Coulomb block¬ 
ade emerge provided that the coupling of the bridge to 
electrodes is sufficiently weak. In this case, there are 
two possible transport channels corresponding to these 
peaks. Omitting from consideration the renormalization 
of the bridge energy Eq due to the temperature gradient 
one may conclude that the best conditions for thermally 
excited electron tunneling in an unbiased system occur 
when Eq = —U and Eq — 0. The curves presented in 
the right panel of the Fig. 2 confirm this conclusion. The 
greatest values of Vth are reached at Eq close to —U 
and/or zero. Slight discrepancies between the predicted 
above and computed values of of Eq providing the best 
conditions for thermally excited transport originate from 
the fact that Eq is renormalized by A9 and this was 
taken into account while computing the results shown in 
the Fig. 2. 

To further analyze the effect of Coulomb interac¬ 
tions between charge carriers on the bridge we con¬ 
sider Vth{A9) at several different values of F assum¬ 
ing that the system is kept at a relatively high tem¬ 
perature {9r > 9k). As presented in the Fig. 3. (left 
panel), in a strongly coupled system where the Coulomb 
blockade peaks are suppressed, Vth takes on small val¬ 
ues. This agrees with an almost featureless character 
of DOS shown in the Fig. 1. However, as the cou¬ 


pling between the QD/molecule and the electrodes weak¬ 
ens, the Coulomb blockade peaks appear. This brings a 
significant increase of thermovoltage values, and makes 
nonmonotonic character of Vth{A9) more distinct. The 
strength of Coloumb interactions also affects the ther¬ 
movoltage behavior. When the charging energy U takes 
on certain values, the dependence of Vth of Ad remains 
nearly linear even when Ad reaches value of the order 
of 9hi, and it may lack a minimum/maximum. An ex¬ 
ample of such a curve is displayed in the right panel of 
the Fig. 3. At other values of U, Vth{A9) becomes an 
nonmonotonic function. The curves shown in the figure 
are plotted assuming that Eq <0. Accordingly, Vth{A9) 
reveals a minimum. The minimum position depends on 
the value of charging energy. The thermovoltage behav¬ 
ior similar to that presented in the right panel of Fig. 3 
occurs when one considers Vth{A9) at different positions 
of the dot energy level provided that the charging energy 
is fixed 351. As discussed in that work, this behavior 


originates from nonmonotonic dependencies of the occu¬ 
pation numbers on the dot of the renormalized energy 
Eq. Actually, {no) is determined by the relationship 


between Eq and U. Therefore, one may observe similar 
features in Vth{A9) behavior either by varying Eq at a 
fixed U or by varying U and keeping Eq fixed as in the 
present work. 

It follows from the previous analysis, that the thermo¬ 
voltage, should strongly depend of the position of the dot 
energy level. As shown in Figs. 4-6, Vth rather sharply 
varies over the interval whose endpoints are located near 
Eq = —U and Eq = 0, respectively. Considering a 
system within the Kondo regime, one sees that as the 
dot level approaches the chemical potentials of electrodes 
from above, Vth reveals a sharp double peak shown in the 
left panel of the Fig. 5. Presumably, this feature appears 
due to the thermally excited hole transport through the 
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FIG. 4: (Color online) The thermovoltage Vth versus Eq and 
temperature differential A9 plotted for a weakly coupled QD 
(r = O.SmeV, U = QmeV) above the Kondo temperature 
\k6R = 0.6mey). 


system which starts when the renormalized dot energy 
level moves sufficiently close to the electrodes chemical 
potentials. The height of the peak increases as the tem¬ 
perature differential A0 enhances. The shape of this 
feature may be treated as that of a single peak with a 
sharp dip imbedded in the middle of the latter. The 
thermovoltage takes on values close to zero on the dip 
bottom. This fall of Vth probably occurs due to the ef¬ 
fect of the Kondo maximum in the electron DOS. The 
presence of the latter gives rise to electron flux which 
counterbalances the hole flux within the corresponding 
range of Eq values. Another and less prominent feature 
appears when the renormalized energy on the dot Eq 
becomes close to —U. It indicates the thermally excited 
electron transport through the channel E = Eq + U. 

Within the Coulomb blockade regime at low tem¬ 
peratures the thermovoltage displays sharp features at 
Eq = —U and Eq = 0 associated with two possible 
channels for thermally excited transport. Both features 
are characterized with derivative-like lineshapes. For the 
considered model representing the QD/molecule by a sin¬ 
gle orbital, the thermovoltage changes sign four times 
over the interval —U<Eo<0, each change indicating 
the change of the kind of charge carriers (electrons/holes) 
predominating in the transport. At higher temperatures 
sharp derivative-like features appearing in the VttiEo) 
lineshape become smoother and broader, as shown in the 
right panel of the Fig. 5. Nevertheless they retain es¬ 
sential characteristics so long as the coupling of the elec¬ 
trodes to the linking QD/molecule remains weak. For 
stronger coupled systems the electron DOS becomes fea¬ 



FIG. 5: (Golor online) Themovoltage versus Eq below the 
Kondo temperature. Gurves are plotted assuming kOn = 
O.lGmeF, U = 6meV for F = 3meV (left panel) and F = 
O.SmeV (right panel). 



FIG. 6: (Color online) Thermovoltage versus Eq above the 
Kondo temperature. Gurves are plotted assuming kOa = 
0.6meF, kOL ~ 0.84meF; F = O.SmeF (left panel) and U = 
QmeV (right panel). 


tureless, and the peaks in the Vth disappear, as shown 
in the Fig. 6. 


IV. THERMOCURRENT 


As mentioned above, properties of thermovoltage in 
nanoscale systems were studied in numerous works. Less 
attention was paid to studies of thermocurrent in spite 
of the fact that Ith is more convenient for measuring 
and modeling [2^. In this Section, we present some re¬ 
sults of theoretical analysis of the thermocurrent behav¬ 
ior in a system consisting of a couple of electrodes linked 
by a QD/molecule which is simulated by a single energy 
level. As before, we assume that the linker is symmet¬ 
rically coupled to the electrodes, and the right electrode 
is cooler than the left one. The thermocurrent is defined 
as follows: 

Ith = nv Or, A0) - nv Or, Ae = 0) (14) 
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FIG. 7: (Color online) Thermocurrent Ith flowing through 
an unbiased QD as a function of the energy Eq and tem¬ 
perature differential A9. The surface is plotted assuming 
r = O.SmeV, kOn = Q.&meV, U = GmeV. 


where V is the bias voltage applied across the system. 
According to this definition, Ith represents the contri¬ 
bution to the tunnel current which originates from ther¬ 
mally excited transport of charge carriers. Employing 
Eq. (IT^ . Ith may be presented in the form: 

Ith ~ J {T{E,9R,Ae)f^{E,eL) 

- T{E, Or, AO = 0)/^ (E, 9l) - ATf^{E, 0^) }dE. 

(15) 

Here, AT is the difference between the electron trans¬ 
mission function computed at nonzero value of the tem¬ 
perature differential A9 and that one computed assum¬ 
ing A9 = 0 : 


AT = T{E,9r,A9)-T{E,9r,A9 = 0). (16) 

In these expressions ((13, dUl), the electron transmission 
is given by Eq. dn. 

The value of Ith is determined by the magnitude and 
polarity of the bias voltage V, and by the temperature 
differential A9. Other factors strongly affecting Eh in¬ 
clude the position of the energy level on the QD/molecule 
Eq, the coupling parameter E, the value of charging en¬ 
ergy U and the average temperature 9. We first consider 
the case of an unbiased system (E = 0). In this case. 
Eh solely originates from the thermally excited flow of 
charge carriers. The thermocurrent dependence of Eq 
and A9 is shown in the Fig. 7. The flow occurs when 
the energy level on the dot is shifted to a position where 




Eg (meV) Eg (meV) 


FIG. 8: (Golor online) Thermocurrent versus Eq for an un¬ 
biased QD within the Goulomb blockade regime below (left 
panel; kdn = 0.16meE) and above (right panel; kOn = 
O.GmeE) Kondo temperature. The curves are plotted as¬ 
suming r = O.SmeV, U = GmeV. 


the renormalized energy Eq (or Eq + U) is close to the 
chemical potential of electrodes. The intensity of the 
flow determines Eh value. It increases as the tempera¬ 
ture differential A9 enhances. This is illustrated in the 
Fig. 8 where the behavior of Eh as a function of Eq is 
shown for several values of A9. The results presented in 
this figure are corresponding to a weakly coupled system 
(r <C U) considered at low {9 < 9k) and moderately 
high [9 > 9k) temperatures. 

At low temperatures, Eh takes on zero value every¬ 
where except close vicinities of two points indicating the 
opening of transport channels for the thermally excited 
transport. In these vicinities, the thermocurrent be¬ 
havior is characterized by sharp and narrow dips whose 
heights increase as A9 enhances. Negative values of Eh 
show that the charge carriers involved in the transport 
process are electrons. At higher temperatures. Eh be¬ 
havior is somewhat modified. As before, the thermocur¬ 
rent takes on nonzero values near Eq = —U and Eq = 0 
which correspond to two transport channels typical for a 
system with a single-state bridge within the Coulomb 
blockade regime. However, the relevant features now 
have derivative-like lineshapes. So, in this case charge 
carriers of both kinds participate in transport via the 
transport channels. At certain values of Eq the two 
fluxes counterbalance each other, so Eh becomes zero. 
Again, the enhancement of temperature differential re¬ 
sults in the increase of the features heights. 

The dependencies of the thermally excited current 
flowing through an unbiased QD/molecule of the tem¬ 
perature differential are strongly influenced by the po¬ 
sition of the energy level Eq. This is illustrated by the 
curves displayed in the left panel of the Fig. 9. These 
curves are plotted for a system within the Coulomb block¬ 
ade regime at a moderately high temperature {9 > 9k). 
Thermally excited transport in such systems was recently 
studied [Sg, and the reported results for Eh qualita- 
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FIG. 9: (Color online) Thermocnrrent as a function of A6' 
within the Colomb blockade regime above Kondo tempera¬ 
ture. The curves are plotted assuming kOR = 0.6meV, T = 
O.SmeV, U = 6meV for an unbiased (left panel) and biased 
(right panel) quantum dot. 


FIG. 10: (Color online) Thermocnrrent as a function of the 
bias voltage in a single-level quantum dot within the 
Kondo regime (left panel) and within the Coulomb block¬ 
ade regime (right panel). Curves are plotted assuming Eo = 
—2meV, U = 6meV for F = 3meV, kOu = O.lGmeV (left 
panel) and for F = O.SmeV, kOu = 0.6meV (right panel). 


tively agree with those presented here. As shown in the 
Fig. 9, Ithi^d) reaches its minimum at a certain value 
of A6 for each considered value of Eq. Besides, some 
of these functions show maxima at small values of A9 
with the subsequent change of sign. As well as in the 
previously discussed case of thermovoltage, the variety 
of Ith{A9) lineshapes originates from the specific non¬ 
monotonic dependencies of the occupation numbers on 
the QD/molecule energy level. 

When a bias voltage is applied across the system, the 
thermocurrent is simultaneously driven by this voltage 
and the temperature differential. The combined effect 
of these two factors determines the Ith behavior. In 
general, A9 remains predominating at small values of 
the bias voltage V. However, for stronger bias, its ef¬ 
fect strengthens to such extent that approaches zero 
regardless of A9 value. In a biased system, the chem¬ 
ical potentials of the electrodes are shifted away from 
previously occupied positions, and the conduction win¬ 
dow opens up. Assuming, as before, that in the un¬ 
biased system chemical potentials associated with both 
electrodes equal zero, we may present them in the form: 
fiL = a|e|y; = —(1 — a)|e|y. In these expressions, 
a is the division coefficient which shows how the bias 
voltage is distributed between the electrodes. In further 
analysis we put a = 1/2 which seems natural for a sym¬ 
metrically coupled system. 

When the conduction window sufficiently broadens to 
contain transport channels, charge carriers start to tun¬ 
nel through the system. However, this transport is af¬ 
fected by the influence of the temperature differential. 
Due to the presence of the latter, two terms in the ex¬ 
pression da do not cancel each other, and Ith may take 
on nonzero values. Curves shown in the right panel of 
Fig. 9 are plotted for several close values of the bias 
voltage. These values are chosen in such a way that the 
energy level on the dot is situated near the boundary of 


the conduction window, either within or barely outside 
it. Therefore, the transport channel at E = Eq appears 
near the shifted chemical potential of the hot electrode. 
This particular position of the dot energy level provides 
the most favorable opportunities for the temperature dif¬ 
ferential A9 to influence the transport by broadening the 
Fermi distribution function associated with the hot elec¬ 
trode. Besides, temperature differential affects the trans¬ 
port by renormalizing the energy Eq. As a result, all 
presented curves display nonlinear lineshapes resembling 
those observed in recent experiments on semiconducting 
quantum dots , and these nonlinear properties of ther¬ 
moelectric transport appear even at small values of A9. 

Some further characteristic properties of thermocur¬ 
rent flowing through biased QD/molecules are illustrated 
in the Fig. 10. In this figure, we show Ith as a func¬ 
tion of the bias voltage for a system within the Kondo 
regime (left panel) as well as within the Coulomb block¬ 
ade regime (right panel). As known, the bias voltage 
application leads to a splitting of the Kondo peak in two 
maxima which are gradually suppressed as the bias in¬ 
creases. The curves displayed in the left panel of the Fig. 
10 are plotted assuming that bias voltage takes on low 
values, and the Kondo peak, although affected, is not 
completely dissolved. Accordingly, each Ith — V curve 
presented in the left panel is characterized by two min¬ 
ima. These features appear due to the presence of two 
peaks in the electron DOS replacing a single Kondo max¬ 
imum in a slightly biased QD. 

When the coupling between the electrodes and the 
linking QD/molecule weakens, and two peaks associated 
with the Coulomb blockade regime emerge, the line- 
shapes of Ith — V characteristics significantly change. 
The interplay between electric and thermal driving forces 
brings several Ith sign reversals, as shown in the right 
panel of Fig. 10. One may suggest that these reversals 

























and accompanying derivative-like features occur as a re¬ 
sult of opening up transport channels at E = Eq, Eq + U 
within the conduction widow. Comparing Figs. 9 and 10, 
one may observe similarities in Ith behavior in unbiased 
and slightly biased systems within the Coulomb block¬ 
ade regime. In both cases derivative like features appear 
in the Ith lineshapes indicating the opening of trans¬ 
port channels, although these features shapes in biased 
systems are distorted due to the bias voltage influence. 
When the bias further increases, electric driving forces 
become predominating. Then the difference between the 
terms in Eq. m diminishes and Ith approaches zero as 
illustrated in Fig. 10. 


V. CONCLUSION 


Finally, we repeat again that thermoelectric proper¬ 
ties of nanoscale systems attract significant interest of 
the research community. The present work was inspired 
by this common interest as well as by recent observa¬ 
tions of strongly nonlinear Seebeck effect in semiconduct¬ 
ing quantum dots bridging two electrodes [i^. Coulomb 
repulsion between electrons on the bridge may signifi¬ 
cantly affect characteristics of thermoelectric transport 
through such systems. The effects of electron-electron 
interactions were theoreticallY a naly zed in several works 
(See e.g. Refs. O [H, [1^ S, [33 • However, in these 
works the analysis was carried out within the Coulomb 
blockade regime (F <C U). Here, we employ a computa¬ 
tional scheme within the EOM approach which enables 
to study thermoelectric transport through QD/molecules 
sufficiently strongly coupled to electrodes. 

We realize that EOM approach used in this work to 
compute relevant Green’s functions has its limitations. 
In principle, better approximations for Green’s functions 
appropriate for studies of electrical and thermal trans¬ 
port in Kondo correlated systems could be derived by 
applying a numerical renormalization group approach. 
Nevertheless, as discussed in an earlier work [20[, EOM 
based approximations for the Green’s functions used in 
the present work bring reasonably good results for trans¬ 
port characteristics assuming that the considered system 
is not too strongly coupled, so that F remains smaller 
that U. 

Applying the suggested computational scheme, we 
showed that the thermovoltage created by temperature 
differential may decrease in magnitude when AO in¬ 
creases. At certain values of A9, Vth even may become 
zero and change its polarity as A9 increases beyond 
these values. This interesting phenomenon was observed 
in experiments and theoretically analyzed for a case of a 
QD within the Goulomb blockade regime [s^. Here, we 
showed that Vth may similarly behave in a strongly cou¬ 
pled system kept below Kondo temperature. However, 
this behavior of thermovoltage is universal in character. 


In general, Vth behavior is determined by interplay of 
several factors including the position of the energy level 
on the considered single-level dot, the average tempera¬ 
ture, the coupling parameter F and the charging energy 
U. We have analyzed some effects of these factors. For 
instance, we showed that at certain values of relevant 
energies, Vth{A9) becomes a monotonic function which 
agrees with results obtained in several earlier works (See 
e.g. Refs. HE]). 

Another important characteristic of thermoelectric 
transport is thermally excited electric current. This 
quantity is especially interesting for it is available for di¬ 
rect measuring in experiments on nanoscale systems. The 
thermoelectric current Ith flowing through an unbiased 
QD/molecule originates from the flux of charge carriers 
from the hot electrode to the cooler one. Within the ac¬ 
cepted model where the dot/molecule is represented by a 
single orbital, there exist two channels for the thermally 
excited tunnel transport. These channels open when ei¬ 
ther the energy on the dot Eq or Eq + U approaches 
the electrodes chemical potentials. Therefore, the value 
of the charging energy U is an important factor affect¬ 
ing Ith behavior. The effect of Coulomb interactions 
strongly depends on the relationship between the energies 
U and F and on the average thermal energy kO. Varying 
these energies, we can move the considered system from 
Kondo regime to Coulomb blockade regime and study 
specific manifestations of Coulomb interactions within 
these regimes. The corresponding analysis is carried out 
in the present work. In a biased QD/molecule, the ther¬ 
mocurrent behavior is determined by the interplay of the 
bias voltage and thermal driving forces. Coulomb inter¬ 
actions affect Ith behavior in biased systems, as well. 
We have analyzed the effect of Coulomb interactions on 
the thermocurrent flowing through a biased QD within 
the accepted model for the latter. We believe that pre¬ 
sented computational method and obtained results may 
be helpful for further studies of thermoelectric properties 
of nanoscale systems. 
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